ON EGYPTIAN FRACTIONS 



MANUEL BELLO HERNANDEZ , MANUEL BENITO , AND EMILIO FERNANDEZ 



Abstract. We find a polynomial in three variables whose values at nonnega- 
tive integers satisfy the Erdos-Straus Conjecture. Although the perfect squares 
are not covered by these values, it allows us to prove that there are arbitrarily 
long sequence of consecutive numbers satisfying the Erdos-Straus Conjecture. 
We conjecture that the values of this polynomial include all the prime numbers 
of the form 4g+5, which is checked up to 10 14 . A greedy-type algorithm to find 
an Erdos-Straus decomposition is also given; the convergence of this algorithm 
is proved for a wide class of numbers. Combining this algorithm with the men- 
tioned polynomial we verify that all the natural numbers n, 2 < n < 2 X 10 14 , 
satisfy the Edos-Straus Conjecture. 



1. Introduction 

One of the most famous conjecture on Egyptian fractions is the Erdos-Straus 
Conjecture (ESC): Given a positive integer n > 2 there exist positive integers 
{x,y,z) £ Z5, such that 

4111 

(1) - = - + - + -. 

n x y z 

In this case we say that n is an Erdos-Straus' number and refer to (JlJ as an Erdos- 
Straus decomposition of — . Sierpihski and Schinzel extent to a more general conje- 
ture replacing 4 in (fTJ) by other fixed positive integer m > 4. The roots of these lie 
in the minimum number of Egyptian fractions needed to decompose a fraction as 
sum of Egyptian fractions (see, for example, [2], [4], pj], [13], [16] and [H]). ESC 
has been verified for all integer up to a bound by many authors: Straus, Bernstein 
P], Yamamoto [19 , Swett [17], etc. Swett has checked ESC for all n < 10 14 . If n is 
an Erdos-Straus' number, then it also holds true for integers which are divisible by 
n; moreover, -. 4 , „ = -4r + ? rrr7i — r^r and if we have the factorization n — a ■ b ■ c, 

' ' 4q+3 q+1 (<3+l)(4<J+3) ' 

! — 1 1 Therefore, ESC should be only checked for primes n = 4q+l. 



a(a-\-b)c b{a+b)c ' 

Because 



(2) 



1 1 1 



abc-a~b a^- a(ac -1)^ ( ac -l)^in 



when be = 1 (mod 4), the numbers n = abc — a — b are Erdos-Straus' numbers. 
Moreover, if we replace the condition be = 1 (mod 4) by be = 1 (mod m), then the 
Sierpihski and Schinzel conjecture holds for n = abc — a — b. Since we are interested 



2000 Mathematics Subject Classification. Primary 11D68; Secondary 11A07,11N32. 
Key words and phrases. Egyptian fraction, Erdos-Straus Conjecture. 

This research was supported in part from 'Ministerio de Ciencia y Tecnologfa', Project 
MTM2009-14668-C02-02. 



2 



M. BELLO-HERNANDEZ, M. BENITO, AND E. FERNANDEZ 



in n = Aq + 1, let n = p(a, /3, 7) be the polynomial p : Z> — > Z given by 

(3) p(a, (3, 7) = (a + l)(4/3 + 3)(4 7 + 3) - (a + 1) - (A/3 + 3). 
This parametric solution of ESC let us to prove the following result: 

Theorem 1. There are arbitrarily long sequence of consecutive numbers satisfying 
Erdds-Straus Conjecture. 

We have checked that the set 

TVi = {n G Z >0 : 3(a,/3,7) e Z> , n = p(a,/3,j)} 

contains all the prime numbers of the form n = Aq + 5, q G Z>o, n < 10 14 . We 
prove, see Lemma [7] below, that Afi does not contain the perfect squares. 

In Section [2] we list several relations for ESC. One interesting conclusion in this 
section is given in Lemma[3l it is impossible to generate "naturally" a finite number 
of congruence classes satisfying ESC which contain all the prime. Section|3]contains 
a proof of Theorem [TJ In Section 0] we prove that the set A/"i does not contain the 
perfect squares and several other additional remarks on ESC. We describe a fast 
algorithm to construct a decomposition of — as sum of three Egyptian fractions 
and we prove the convergence of that algorithm for a large class of numbers in the 
last Section [5] 

2. Auxiliary relations 

The Erdos-Straus Conjecture shows an arithmetic property of the natural num- 
bers (see Lemma E] below). We can split any fraction ^ as sum of two Egyptian 
fraction but the same does not hold for — for all n G Z>o- It is very well known 
that the equation ^ = i + ^ has solutions (x, y) G Z?, if and only if n has a divisor 
m with either m = or m = 2 (mod 3); it is equivalent to n ^ 1 (mod 6). 

The following Lemma is very well known; see, for example [13], p. 287, or [19] . 
We give a proof for an easy reading. 

Lemma 1. Let n be a prime; n is an Erdds-Straus' number if and only if there 
exists (a, b, c, d) G Z 4 such that some of the following conditions holds: 

(4) (Aabc - l)d = (a + b)n, 

(5) (Aabc - l)d = an + b. 

Proof. If (Aabc — l)d = (a + b)n or (Aabc — l)d = an + b holds, dividing these 
equations by abedn, we have respectively 

(6) 4 - 1 + 1 + 1 

n aben bed aaV 

(7) - - — ^ + — + 1 

n aben bed acdn 

On the other hand, if ([1]) holds, then Axyz = n(xy + yz + zx). Since n is prime, n 
divides x, y or z. Of course, n does not divide all the three numbers x, y, z because 
this trivially yields the contradiction 4 = 77^ + ^7^ + 77^ with x/n, y/n, z/n positive 
integers. Hence, we can assume without lost of generality x — an. Thus (p} is 
equivalent to 

4a - 1 1 1 1 1 1 

( 8 ) = 7. + -^ 



y z na (4a — l)y (4a — l)z 
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Since n is prime and (4a — l,a) = 1, we have two cases: (4a — l,n) — 1, and 
(4a — 1, n) = n. 

In the first case, there exist 01,02,03 £ Z + such that a = 010203, (nai,ci2) = 1, 
and 

(9) (4a — l)y = na\{na\ + 02)03, (4a — l)z — 02(1101 + 02)03. 

Because (no, 4a — 1) = (71010203, 4a — 1) = 1, there exist a, ft such that y = 
anaid3, z — ,80203. From (j9]), we have (4a— 1)q = (7101 + 02) = (4a— l)ft, it means 
a = ft. Therefore, A — ai, B — 02, C — 03, and D = a satisfy 

(AABC - l)D = {nA + B), 

Now we consider the second case (4a — 1, n) — n. So, there exists j such that 

(10) 4o - 1 = jn. 

Setting this expression in ((5]), we have - = J- + -i. Thus, there exist oi, 02, 03 G Z + 
such that a = 010203, (01,02) = 1 and 

(11) jy = ai(ai + a 2 )a 3 , jz = a 2 (ai + a 2 )a 3 . 

Since (j, a) = (j, 010203) = 1, there exist a, ft such that y = 00103, z = fta^a-z- 
Setting these expression in (TTTI) . we have ja = a\ +02 = jft, it is a = ft. Multiplying 
in (|10[) by a, we obtain A = a\, B = 02, C — 03, D — a such that 

(AABC - l)D = [A + B)n. 

□ 

The relations ((4]) and ([5]) are equivalent to © and ([7]), respectively. Since in the 
first decomposition of ^, n divides one denominator but is coprime to the others, 
while in the second it is only coprime to one of them, these are respectively refereed 
as Type I and Type II decomposition of — (see [2]). Observe that it is not possible 
for all three of x, y, z in (fTJ) to be divisible by n. 

Changing variables in (|4]) and ([5]) the following lemma follows. 

Lemma 2. (1) Let n be a prime. There exist positive integers a, 6, c, d such 
that holds if and only if there are positive integers a, ft, 7, S such that 

(12) Sn = (4aft-fS - 1) - 4a 2 7 . 

(2) Let n 6 Z, n > 2. There exist positive integers a,b,c,d such that {5[) holds 
if and only if there are positive integers a,ft,j,6 satisfying 

(13) n = (4a/3 7 - 1)5 - 4ft 2 j. 

Proof. (1) If n is a prime and (0| holds, then d divides a + b. Set e = 

hence, b = de — a and (@| yields en = (Aacde — 1) — 4a 2 c. Setting d = e, 
a = a, ft = d, and 7 = c, we obtain (1121) . 
(2) The relation $5$ is equivalent to b + d is divisible by a and n + s = 46cd, 
where s = <4> d = as — 6. Setting a = a, /3 = 6, 7 = c, 8 = s we obtain 
immediately fp~3|) . 

□ 

Equation (|12p implies that on is an Erdos-Straus' number with Type I and Type 
II decomposition, but it is trivial because if n = 1 modulo 4 it follows S = — 1 modulo 
4. The parametric relation (|13j) appears implicitly in |18j . Setting ft = 7 = 1 in 
(fT3|) . it follows that if n + 4 has a divisor congruent to 3 modulo 4, then n is an 
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Erdos-Straus' number. So by the Landau prime ideal theorem (see [9] or [12] pp. 
266-268) the set of positive integers which are no Erdos-Straus' numbers has zero 
density. A sharp estimate on the density of Erdos-Straus' numbers can be found in 
PI (see also [2], [6], [7], and 0). 

In [19] Yamamoto (see also [15]) proves that a perfect square n does not satisfy 
neither (j4]) nor ([5]) with some restriction in the parameters (see the paragraphs 
at the begging of Subsection 14.21 below), so fixing the parameters a,b,c,d in these 
equations we can not generated a complete residue system. The following lemma 
emphasis in this remark without use that result of Yamamoto. 

Lemma 3. We can not generate a finite number of congruence classes containing 
all the primes of the form 4g + 5 fixing in or (0) three of the four parameters 
in a finite subset ofl? >Vl and the remain parameter free in Z >0 . 

Proof. Let S a — {(6, c, d) G Z?, } denote a finite subset of Z5, of values taken 
by (b,c,d). From (U) or ([5]), when a is free taking any positive number, each 
fixed vector (6, c, d) G S a may generate a residue class of Erdos-Straus' numbers n. 
Equivalently we define Sb,S c and Sd- We fixe four such subset S a , Sb, S c , Sd and 
first prove that there exist infinite prime numbers n which can not be generated by 
([5]) with parameters (a, b, c, d) with three of them in the corresponding set S a , Sb, S Cl 
or Sd- 

The numbers n generate by ([5]) with (b, c, d) G S a are a finite set because a\(b+d). 
Moreover, observe that in the proof of (|13p in Lemma [21 the numbers n generate 
by ([5]) or equivalently by (fl3|) with (b, c, d) G S c are given by 

n = 4(a5 - /?)/3 7 - S = Abed - - + - , (a, b, d) G S c , 

a 

where a = a, = b, 7 = c, and 5 — 2±2, Taking T c = \cm{bd : (a,b,d) G 
S c }, the numbers in {AT c t + 1 : t £ Z>o} (in particular the prime numbers in 
{AT c t + 1 : t G Z >0 }) can not be generated by ([13]) with (a,b,d) € 5 C , here 
lcm{6d : (a, 6, d) € S'c} denotes the least common multiple of the numbers bd when 
(o, b, d) runs through the set S c - Suppose for contradiction that for a given t G Z>o 
there exist (ao, 60, do) G S c and c G Z>o such that 

4T c i + 1 = 4& cd - 5°±* o 1 = 4& d (c - e) - , 

where e G Z >0 , but it is a contradiction since 1 is not an Erdos-Straus' number. 

Using the same arguments and notations as before, when (a, 6, c) G Sd, then 
according to the proof of (fT3"]) in Lemma [5] (a, /5, 7) lies in a finite subset S5 of 
Z> , and (5 G Z>o- Taking = lcm{(4o;/?7 — 1) : (a, j3, 7) G S5}, the numbers in 
{4Tdt + 1 : t G Z >0 } can not be generated by (fl"3]) with (a, 6, c) G S^. Suppose for 
ate Z>o there exists (ao,/3o,7o) G 5s and S G Z>o such that 

AT d t + 1 = (4a /?o7o - 1)<5 - 4/?g 7o ^ 1 = (4a 0/ 3o7o - 1)(* - 4e) - 4^ 7o 

for some e, again it is not possible because 1 is not an Erdos-Straus' number. From 
the symmetry of (| L3[) in b and d the same conclusion holds for Sb for numbers 
{ATbt + 1 : t G Z >0 }, where is defined as T^. Of course all the numbers 

{4T b T c T d t + l:t GZ >0 } 

can not be generated by parameters (a, 6, c, d) with three of them in the correspond- 
ing set S a ,Sb, S c or S d - 
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The same arguments work for (j4|) . For example, one of the parameters a or b in 
(J4j) can not take values in an infinite subset of positive integers while other three 
belong to a finite subset. In fact, observe that Q is symmetric in a and &, and 
according to the proof of (|12l) in Lemma [3J the number e = divides 1 + 4a 2 c, 
so that if (a, c, d) runs over a finite set, the numbers of divisors of 1 + 4a 2 c is finite. 
Moreover, the numbers n generate by ((3]) or equivalently by (TT21 with (b, c, d) £ 5 C 
do not contain {4T^ + 1 : t £ Z>o} where T' c — lcm{aM : (a, b, d) £ S'c}- Suppose 
for contradiction that for a given t there exist (ao,bo,do) £ >Sc and c £ Z>o such 
that 

(a + b ){4T' c t + 1) = (4a 6 c - l)d (<x + 6o)l = (4a 6 (c - e) - l)do 
where e = (ggjj°)Zkj but j g a contradiction since 1 is not an Erdos-Straus' 

a b do ' 

number. For (a, b, c) £ Sd, we can not generate the numbers 4T' d t + 1, where 
= lcm{(4a&e - 1) : (a, b, c) £ Sd}. 

Therefore, equations (U]) and ([5} can not generate all the numbers 
(14) {AT b T c T d T' c T' d t + 1 : t £ Z >0 } 

except possibly a finite number of them with three parameters of a, 6, c, c£ in the 
corresponding set S a , S&, S c or S d . 

□ 

Remark 2.1. Note that using Yamamoto result we can actually prove that equations 
(gl) and J5]) can not generate all the numbers {4T b T c T d T' c T' d t + n : t £ Z >0 } 
except possibly a finite number of them with three parameters of a,b,c,d in the 
corresponding finite set S a ,Sb,S c or Sd, where no is a quadratic residue modulo 
AT b T c T d T' c T' d . 

Many papers devote special attention to Type II decomposition since parametric 
solutions of ESC are easy obtained for this case, see (TT51) . Our following lemma 
allows us to find a parametric solution of Type I decomposition. 

Lemma 4. Let n £ Z>o- There exist positive integers a,b,c,d such that holds 
if and only if there are x, t, A such that 

X n+t r- rjj 
J £ A>0; 



^ £Z 



(15) 
are satisfied. 

Proof. Assume that there exist positive integers x, t, A, y, z such that 
x n + t n + A 

= y, — z x n + t = y X, A = 4 z xt — n 

A y ' Axt y ' 

(x + y)n = (Axy z — l)t. 
Therefore, taking a — x,b — y. c = z, d = t we have (Aabc — l)d = (a + b)n. □ 



Remark 2.2. If n is an Erdos-Straus' number and (|15)) is satisfied for positive 
integers x,t,X, then for all j £ Z>o, A*" = n + ixtXj is also an Erdos-Straus' 
number. In fact, 

f xN+t = x(n+A X tX 3 )+t =£ n±t + ^ g z+ 
1 JV+A _ (n+4xt\ j)+X _ n +\ , x ■ c ^ 
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In particular, the set of Erdos-Straus' numbers is an open set in Furstenberg's 
topology (see [13] . p. 34). 

Remark 2.3. The values 

n = p(a, /3, 7 ) = {a + l)(4/3 + 3)(4 7 + 3) - (a + 1) - (4/3 + 3) 

satisfy HU) with x = 1, t = a + 1, and A = 4/3 + 3. 

The parametric relation ([3]) for p = 4g + 5 is useful to rewrite for q. Since 
p(a, p, 7 ) = 4«(a, ft 7 ) + 5 = 4 (((4/3 + 3) 7 + (3/3 + 2))(a + 1) - {fi + 2)) + 5, 

the following lemma follows: 

Lemma 5. If there exists (a, ft 7) G Z>o smc/i £/ia£ 

(16) q = q(a, ft 7) = ((4/3 + 3)7 + (3/3 + 2))(a + 1) - <fi + 2), 

f/ien the Erdos-Straus Conjecture holds for p = 4g + 5. 

3. Consecutive numbers satisfying ESC 

Next, using the Chinese Remainder Theorem and Lemma[5]we prove Theorem[T] 
Actually, we prove that there exist arbitrarily long sequence of consecutive residues 
class such that ^ has Type I decomposition for all n in this residues class. 

Proof of Theorem [7J Of course, to get Theorem [T] it is enough to consider 
"consecutive" numbers n = 1 mod 4. Let N be an arbitrary positive integer and let 
A be a subset of Z>o containing N consecutive nonnegative integers, for example 
A = {0, 1, . . . , N — 1}. If {/3i, (32} C A, then the greatest common divisor of 4/3i + 3 
and 4/32 + 3, (4/3i + 3, 4/3 2 + 3), is an odd number and 

(4/3! + 3,4/3 2 + 3)1(4/3! + 3 - (4/3 2 + 3)) = 4(/3i - /3 2 ) 

=>• (4A + 3, 4/3 2 + 3)|3(/3i - /3 2 ) = 3ft + 2 - (3/3 2 + 2). 

Thus, by the Chinese Remainder Theorem, there is a natural number T such that 

T = 3/3 j + 2 (mod 4ft- + 3), Vft G A, 

i.e., there exist positive integers 7^ such that 

T = {4(3, + 3) 7j + 3ft + 2, Vft G A. 

According to Lemma [S] all n = Aq + 5 with g in the consecutive residue class 
— (ft + 2) (mod T), ft G A, satisfy the Erdos-Straus conjecture. 

□ 

Remark 3.1. Moreover, we can also prove that there exists arbitrarily long sequence 
of consecutive numbers n such that ^ has Type II decomposition. Let A be a set 
of natural numbers containing consecutive numbers, for each a G A, we choose 
natural numbers /3(a) and 7(a) such that a = /3(o) 2 7(a) (such representation is 
unique taking 7(a) free of squares). We chose T as the least common multiple of 
(4/3(a)7(a) — 1), where a takes the values in A, i.e. 

T = lcm{(4fta)7(a)-l):aGA}, and . f ^ ^ £ 

then, using (fl3|) in Lemma [21 it follows that all the fractions T<5 ^ 4a , a G A, have 
Type II decomposition. 
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4. Additional remarks 

4.1. (/ Conjecture. 

Lemma 6. ESC is true if and only if for each q € Z >0 there exists (x, y, z) £ 
such that one of the following relations holds 

(17) q = 1 + 3x + 3y + 4xy, 

(18) q = 5 + 5.x + 5y + 4xy, 

(19) q = q(x,y,z). 
where q{x,y,z) is defined by &16]) . 

Proof. Because 

q = 1 + 3x + 3y + 4xy & 4g + 5 = (4.x + 3)(4y + 3), 

q = 5 + 5x + 5y + 4xy <s> Aq + 5 = (4(x + 1) + l)(4(y + 1) + 1), 
all q for composite numbers Aq+5 satisfy (fTT)) or (|T5| . We also have for q = q(x, y, z) 
define by HI]), 

p(x,y,z)-5 
g = g(a;, y, z) = . 

Hence if all q € Z>o satisfies (fT7)) . (fT8| or (fT9]) , the value of g for all prime numbers 
of the form 4g + 5 lies in M\ and the proof is concluded using Lemma [SJ □ 

Remark 4.1 . We have verified that N\ contains all primes of the form n = 4q + 5 < 
10 14 . To this aim we have generated the classes of equivalence of numbers contained 
in A/i translating t units each variable; i.e. 

N x = p(x, y, z) + h(y, z)t ^ p(x, y, z) + 4 (% + 3)( ^ + 3) ~ \ 

N y = p(x, y, z) + f 2 (x, z)t d = p(x, y, z) + 4 ((as + l)(4z + 3) - l)t, 

N z - p(a;, y, z) + f 3 (x, y)t = p(x, y, z) + 4{x + l)(4y + 3)t. 

Starting with x,y,z E {0, 1} we obtain the following equivalence class 

5 + 8t, 5 + 12t, 13 + 20t, 17 + 20t, 13 + 28t, 37 + 52t, f e Z> , 

then we sieve prime numbers congruent to 1 modulo 4 in these congruence classes. 
We are checked by computer calculations that the remaining primes < 10 14 belong 

to Mi. 

Conjecture. We conjecture that all the natural number q € Z>o can be writ- 
ten as one of the three above relations (|17p . (fT8|) or (|19p . We refer it as the 
g— Conjecture. According to the lemma above it is equivalent to all prime numbers 
of the form 4g + 5 lie in Mi ■ 

4.2. Mi does not contain perfect squares. Here we prove that Mi does not 
include perfect squares, we need to use Jacobi's symbol. Let p be an odd prime 
and n£Z with (n,p) = 1, Legendre's symbol is defined by 

' n \ J 1, if n is a quadratic residue mod p, 
pi 1—1, if n is a quadratic non- residue mod p. 

Let the standard factorization of m be p\p2 ■ ■ -pk where the p r may be repeated. 
If (n, m) — 1 then Jacobi's symbol is defined in term of the Legendre symbol by 
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(m) = rij=i fpO ■ ^ e following properties of the Jacobi symbol are very well 
known: Let m and m! be positive odd integers, (i.) If n = n' (mod m) and 
(n,m) = 1, then (£) = (ii.) If (n,m) = (n,m<) = 1, then (£) (£) = 

(^)- (Hi-) K( = (n» = 1, then (£) (£) = (aal). (iv.) If „ = 3 (mod 

4), (— 1/n) = — 1. (v.) Law of reciprocity. Let m and n be odd coprime. Then 
(m) (^) = ( — l) 1 ^ --12- . (vi.) If 7i, d are coprime positive integers and m satisfies 
n = — m (mod d), then (^) = f^-) , see [10], p. 305. The Kronecker symbol is 
defined in term of Jacobi's symbol. It satisfies (^) = for (n, m) > 2. 

In |19j . Yamamoto observes that numbers n, satisfying (4et6c — l)d = (a + 6)n 
for some positive integers a, 6, c, d with (n, aod) = 1, are no perfect square. Observe 
the symmetry of (4abc — l)d = (a + b)n in a and b and Yamamoto proves, using our 
notations, (jh) = —1 for the Kronecker symbol. Moreover, if (4aoc—l)d = (a+b)n, 
then (4aoc — l)d' = (a + b)n 2 with d' = dn. Our class A/i contains numbers n which 
are not included in Yamamoto class. For example, it is the case of n = 2009, where 
a = 1, b = 293, c = 12, and d = 42, meanwhile a + 1 = 42, 4/3 + 3 = 7, and 
47 + 3 = 7. In fact, we have 

2009 = 42-7-7 - 42 - 7, 
2009(1 + 293) = (4 • 1 • 293 - 1)42, 
(2009,293-42) = 7. 

Lemma 7. A/i does not contain perfect squares. 

Proof. Let n £ A/i. Then there exist nonnegative integers a, /3, 7 such that 

n + (4/3 + 3) = (a + 1) ((4/3 + 3)(4 7 + 3) - 1) = f {a + l)r. 

where r = (4/3 + 3) (47 + 3) — 1. Because of (4/3 + 3, r) = 1, we also have (n, r) = 1 
and (n + t, t) = 1. Using reciprocity law for Jacobi's symbol, we obtain 

/ n \ /n + r\ , , n -i n+x-i 

(20) krr- =(-!)^ — a ~ = i, 



here we have used n = 1 (mod 4). Since n + r = r (mod n), we have f 1 ^ 2 -) = (£) ■ 
Taking into account n = —(4/3 + 3) (mod r) and property (vi) for the Jacobi's 

symbol cited above, we have (~) = ( 4 p +3 J ■ Since t = — 1 (mod 4/3 + 3) and using 
Property (iv), we obtain we have 

-^7) = —1, it yields n is not a perfect square. □ 



Since 4(n 2 + n — 1) + 5 = (2n + l) 2 , and 4g(a, /3, 7) + 5 = /3, 7), where 
p(a, /?, 7) and q(a,/3,j) are given by ([3]) and (fl6|) respectively, the Lemma [7] gives 
immediately the following result: 

Corollary 2. The numbers n 2 +n + /3 + 1, (n,/3) € ^> : /lave no divisors congruent 
to 3/3 + 2 modulo 4/3 + 3. 

4.3. g— Strong Conjecture. Lemma [7] and computer calculations let us to for- 
mulate the following conjecture which implies ESC and has been checked up to 
2 x 10 12 . 
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q— Strong Conjecture: The positive integers Z>o is disjoint union of the three 
sets: 



• A = {q : 3n 6 Z>o, q = n 2 + n — 1}. 



£? = {q : 3(a, /3,7) <E Z2, , q = <j(a, ,8,7)}, where q(a,P, r y) is given by ([16 



C = {25, 115, 145, 199, 659, 731, 739, 925, 1195, 1235, 2381, 3259, 3365, 

3709, 4705, 6325, 8989, 15095, 27991, 39239, 62129, 174641, 279199, 

281735, 310771, 404629, 1308259, 1822105, 2083075}. 

Of course this conjecture implies the q— Conjecture. 

5. A Greedy-Type Algorithm for the Erdos-Straus Decomposition 

Let us describe an algorithm to construct a decomposition of — as sum of at 
most three Egyptian fractions. Using this algorithm we have checked ESC for n up 
to 2 x 10 14 . In Lemma [5] bellow we prove that the following algorithm converges 
for all numbers n = abc — a — b with b = c mod 4. Meantime in Lemma [H] we give 
some conditions which characterize the convergence of this algorithm. 

A Greedy- Type Algorithm for the Erdos-Straus Decomposition: given 

jieN,n>2, 

Step 1. Set j = 1 and q = |_f J , the integer part of j. 
Step 2. Set Xj=q + j, Kj ■ = £ - ±, and yj = Hi-]. 

Step 3. If ^ — 7- — ^7" — Ois true, the algorithm stops else we ask if Zj — ±_J___i_ £ 



V] 

1 is true, then the algorithm stops also and - = — + — + — else we set 

n Xj yj zj 

j i ■H j ' + 1 and return to Step 2. 

Here y = [~— ] denotes the ceiling function evaluated at — . 

We call this algorithm "The Greedy-type Algorithm for Erdos-Straus decom- 
position." Observe that the algorithm finds the greedy decomposition of , 
j = 1,2,..., and also stops when the decomposition is the sum of two Egyptian 
fractions. In fact, if Tj = —nxj (mod 4j — 1) with fj G [1, 4j — 2], then 

4 j - 1 1 Vj 
nxj yj nxjyj 

The algorithm stops when rj divides nxjyj and zj = nx £ Vi . 

Lemma 8. The greedy-type algorithm converges for all n = abc — a — b, with be = 1 
mod 4. 

Proof. If b = c = 1 mod 4, then n = abc — a — b = 3 mod 4, in this case ^ has 
a decomposition as sum of two Egyptian fraction. So we only have to consider 
b = c = 3 mod 4. 

First, we observe that if (x, z) £ Z?, , then 1/x is the Egyptian fraction nearest 
to — + - and less than this number if and only if 

X z J 

(21) x(x-l)<z. 
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Moreover, we have 

4411 1 
(22) — = = l l 

V ' n abc-a-b a^- o(oc-l)^l ( ac -l)^n' 

and 

be — 1 . .be — 1 . .be — 1 
a — - — < a(ac — 1) — - — < (ac — 1) — - — p 

4 v 4 4 

since c = Az + 3 > 3, ac — 1 > 2a, and n = abc — a — b > a b(ac — 1) > 2a. 
According to (|2ip and (|2"2"j) . it is enough to prove 

be - 1 / be - 1 \ . , 6c - 1 

a — - — a — ; 1 < (ac — 1) — - — n 

4 V 4 / ~ 4 

a ^a — 1^ < (ac — l)(a6c — a — 6). 

Observe that both relations a (a^ 1 ^ — l) = 6a 2 c/4 — a 2 /4 — a and (ac — 1) 2 6 — 
a(ac — 1) are linear function in 6, so comparing their leading coefficients and their 
values at b = 1, the above inequality follows. 

□ 

Remark 5.1. We are checked that the algorithm given above also works for Sicrpihski's 
Conjecture, ^ = i + i + i, taking g = [jj and making other few trivial changes, 
for n up to 10 10 . 

If n = 4q + 1 and (4g + l)(g + j) = (s + l)(4j - 1) + r - (4j - 1), then 
4 1 1 4j - 1 - r 



(23) 

4g + l q + j s + 1 (4g + l)(g + j)(s + l) 

and 

(24) (4q+l) 2 (q + j) 2 = (4g+l)(g+i)(s + l)(4j-l) + (r-(4i-l))(4g+l)(g+i). 
So the first part of following result follows: 

Lemma 9. Set n = Aq + 1 G Z>o- £e£ s = s(q,j) and r = r(q,j) be such that 
(Aq + l)(q + j) = s(Aj - 1) + r, < r < 4j - 2 for some j £ Z >0 . T/ien ifte 
greedy-type algorithm converges for n = Aq + 1 j/ and ora/y i/ one of t/ie following 
equivalent statement holds: 

(1) (4j -l)-r divides to (Aq + l)(g + j)(s + 1). 

(2) (4j - 1) -r divides to {Aq + l) 2 (g + j) 2 . 

Proof. According to (|23|) . it is enough to check that if (4j — 1) — r divides to 
(4g + l) 2 (g + j) 2 , then (4j - 1) - r also divides to (4g + l)(g + j)(s + 1). Let 
5 = (Aj — 1, s + 1), so we have (a, /3) = 1, where a = J 7 and (3 = S ^ L - Since J is 
odd, (4g + l)(g + j) = (4g + l) 4g+1+ 4 (4j - 1} , and 
(25) 

(4g+l)(g+j)(s+l) = ( S +l) 2 (4j-l)-(4/c-l-r)(s+l) = {s+l) 2 Sa-{s+l)5(a-P), 

if 5 divides (4g+ l)(g + j)(s+ 1), then S 2 divides (4g+ l)(q + j)(s + 1). On the other 
hand, using (|24|) and the hypothesis we obtain a — f3 divides (4g + l)(g + j)(s + 1). 
Therefore, from (|25p we conclude that (4j — 1) — r = 6(a — /3) divides (4g + l)(g + 
j)(s + 1), and the lemma is proved. □ 
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Remark 5.2. For theoretical considerations is useful to observe that according to our 
notations (4q+l)(q+j) = s(4j-l)+r <S> (4<7+l) 2 + (4<7+l)(4j-l) = 4s(4j-l)+4r, 

so we have (4q + l) 2 = (4s - (4q + l))(4j - 1) + 4r =' (4k - l)(4j - 1) + 4r and 
(4q + l) 2 = 4r mod 4j — 1. 

Remark 5.3. From statement 2 in the lemma above it follows that if n = 4q + 1 is 
a prime, the greedy-type algorithm stops at a step j with 4j — 1 < 4q + 1 if and 
only if 4j — 1 — r divides (q + j) 2 . 

Remark 5.4. Given any j G Z>o there exists a prime number n = rij such that the 
number of steps in the algorithm above is larger than j. In fact, if 

n d = 41cm{3, 7. . . ,4j - 1, 2, 5, . . . , 3j - l}t + 1 =' 4 qj + 1 

is a prime for a given i G Z >0 , then for each /c < j, we have 

(iqj + l)(qj +k) = k, mod 4fc - 1, 

and since (k, 3fc — 1) = 1, 

4fe - 1 - k = 3k - 1 j (Aqj + 1) 2 (^- + kf . 

Therefore, according to the Lemma[9l the greedy-type algorithm does not converge 
in the first j steps for rij. 

6. Appendix. Computer Programs 

In this section we include some programs which are used to check the Erdos- 
Straus Conjecture, q— Conjecture and the q~ Strong Conjetures. For an indepen- 
dent reading we set again some facts already cited. 

6.1. A computer program checking the q —Strong Conjecture. Lemma [7| 
and computer calculations let us to formulate the following conjecture which implies 
the ESC and has been checked up to 2 x 10 12 : 

q— Strong Conjecture: The set of positive integers, Z>o, is disjoint union of 
the following three sets: 

• A = {q : 3n 6 Z >0 , q = n 2 + n - 1}. 

• B = {q: 3(a, /3, 7) 6 Z| , q = q(a, (3, 7)}, where 

(26) q = q(a, 0, 7) = ((4/3 + 3) 7 + (3/3 + 2))(a + 1) - (/3 + 2), 

C = {25, 115, 145, 199, 659, 731, 739, 925, 1195, 1235, 2381, 3259, 3365, 

3709, 4705, 6325, 8989, 15095, 27991, 39239, 62129, 174641, 279199, 

281735, 310771, 404629, 1308259, 1822105, 2083075}. 

Of course this conjecture implies the q— Conjecture, see below and LemmalU Next 
we describe our algorithm: 

Step 1: We generate a matrix w of equivalence classes in the set B using the 
function cribata of two arguments. The first column of w are the moduli 
and the second are their corresponding rests. We only use moduli which are 
divisors of a given number. The values in q of two of the three parameters 
a, /3, and 7, or their associated new parameters after a change of variable, 
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are setting up to a bound. We use several equivalent expressions to the 
function (l26l): 



q = Axyz + 3xy + 3xz + 2x + Ayz + 2y + 3z, 

g = x(4yz + 3y + 3z + 2) + 4yz + 2t/ + 3z, 

q = y(Axz + 3x + Az + 2) + 3xz + 2x + 3z, 

q = z(Axy + 3x + Ay + 3) + 3xy + 2x + 2y, 

Aq + 6 + x = (Ay + 3)(4(> + \)z + 3x + 2), 

q + y + 2=(x + l)((4y + 3)z + 3y + 2), 

(3 + 4z)q + 4 + 5z = ((3 + 4z)a; + 4z + 2)((3 + 4z)y + 3z + 2), 

and others which are easy to find after a change of variables, setting y = z+a 
in (|27|). we have 

g = x(4z 2 + z(4a + 6) + 3a + 2) + 4z 2 + z(4a + 5) + 2a, 

doing z = tl ~°~ 1 in (1271) with it, a of different parity, — u < a < u, 



2 

, , 2 a + u + 1 
q = u + u — 1 — a h cc(u + u — a") 

setting a = u — 2d — 1, 

= (it - d- l)(4d + 2) + 3d + x{u{Ad + 3) - (2d + l) 2 ), < d < - 



2 1 
u- 1 



for a = — u + 2d + 1, we have 

q=(u-d- l)(4d + 3) + 2d + x(u(4d + 3) - (2d + l) 2 ), < d < 

doing x = s — t, y = t — z in (|27|) . we obtain 

g = s(-4z 2 + 4<z + 3i + 2) + 4iz 2 - 4z 2 - 4t 2 z + Atz + z - 3t 2 , 
doing the same x — s ~ t, z — t ~ y in (|27l) . 

g = (-4y 2 + 4ty + 3t + 2)s + Aty 2 - 4y 2 - At 2 y + Aty - y - 3t 2 + t, 
with y = s — t, z = t — x in (|27l) . 
g = (-4x 2 +Atx-x + At + 2)s + Atx 2 - 3x 2 - At 2 x + Atx - x - At 2 + t, 
doing y = s — t, x = t — z in J27J), 

q = s(-Az 2 + Atz + z + 3t + 2) + Atz 2 - 3z 2 - At 2 z + 2tz + z - 3t 2 , 
setting z = s — t, y = t — x in (|2"T|) , 

q = s(-Ax 2 + Atx - x + At + 3 + Atx 2 - 3x 2 - At 2 x + Atx - At 2 - i, 
with z = s — t, x = t — y in (|27p . 

g = s(-4y 2 + Aty + y + 3t + 3) + Aty 2 - 3y 2 - At 2 y + 2ty - 3t 2 - t. 

Step 2: Next, we sieve g in the residue classes in the vector w using the 
function cricu. 

The steps 1 and 2 are done two times. First, we use moduli which are 
divisors of the product of all primes up to 19. Next, we add the factor 23 
obtaining 

tc = 2 x 3 x 5 x 7 x 11 x 13 x 17 x 19 x 23. 
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Step 3: As a by-product of steps above we have many rests modulo tc which 
are not sieved, all these are saved in the vector v. The function adjunta 
attaches to a matrix a vector, this is also used in the Step 1 generating 
the matrix w. At the last, we use the function fasf which combines square 
checking with functions cri throughout the equivalent classes in v modulo 
tc. 

Next, we show Sage code of our program, although we write also our algorithm 
in UBASIC and Pari GP. The Sage code is more transparent and easy to read. 
Moreover, the program have been run on the operating systems: MS-DOS, Linux, 
and OS X. 

def adjunta(s,r,u) : 
l=len(u) 
h=0 

for i in xsrange(l): 
if s°/„u [i] [0] ==0 : 
c=gcd(s,u[i] [0]) 
if r%c==u[i] [1] : 
h=l 

if h==0: 
1=1+1 

u. append ( [s ,r] ) 
return (u) 

def cribata(t ,m) : 

u=[[2,0] , [3,0] , [7,0]] 
ml=m+l 

for b in xsrange (ml) : 

for a in xsrange (b+1) : 
s =4* a *b+3*a+3*b+2 
if t°/„s==0: 

r=4*a*b+2*a+3*b 

u=adjunta(s,r,u) 
s=4*a*b+3*a+3*b+2 
if t'/„s==0: 

r=4*a*b+3*a+2*b 

u=adjunta(s ,r ,u) 
s=4*a*b+3*a+4*b+2 
if t7.s==0: 

r=3*a*b+2*a+3*b 

u= ad j unt a ( s , r , u ) 
s=4*a*b+4*a+3*b+2 
if t°/„s==0: 

r=3*a*b+3*a+2*b 

u=adjunta(s,r,u) 
s =4*a*b+3*a+4*b+3 
if t'/„s==0: 

r=3*a*b+2*a+2*b 
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u=ad j tint a ( s , r , u) 
s =4* a *b+4*a+3*b+3 
if t%s==0: 

r=3*a*b+2*a+2*b 

u=ad j unt a ( s , r , u) 
s=4*a*b-4*a~2+3*b+2 
if t%s==0: 

r=4*a~2*b-4*a~2-4*b~2*a+4*a*b-a-3*b~2+b 

u=ad j unt a ( s , r , u) 
s=4*a*b-4*a~2+3*b+2 
if t%s==0: 

r=4*a~2*b-4*a~2-4*b~2*a+4*a*b+a-3*b~2 

u=adjunta(s ,r ,u) 
s=4*a*b-4*a~2+4*b+2-a 
if t%s==0: 

r=4*a~2*b-3*a~2-4*b~2*a+4*a*b-a-4*b~2+b 

u=ad j unt a ( s , r , u) 
s=4*a*b-4*a~2+a+3*b+3 
if t%s==0: 

r=4*a~2*b-3*a~2-4*b~2*a+2*a*b-3*b~2-b 

if r>0: 

u=adjunta(s ,r ,u) 
s=4*a*b-4*a~2+a+3*b+2 
if t%s==0: 

r=4*a~2*b-3*a~2-4*b~2*a+2*a*b-3*b~2+a 
if r>0: 

u=ad j unt a ( s , r , u) 
s=4*a*b-4*a~2-a+4*b+3 
if t°/.s==0: 

r=4*a~2*b-3*a~2-4*b~2*a+4*a*b-4*b~2-b 
if r>0: 

u=adjunta(s ,r ,u) 
s=-4*a~2+4*a*b-4*a+3*b-l 
if s>0: 

if t°/.s==0: 

r=- 4*a~2 + 4*a*b - 5*a + 3*b - 3 
if r>0: 

u=adjunta(s , r ,u) 
s=-4*a~2+4*a*b-4*a+3*b-l 
if s>0: 

if t7.s==0: 

r= - 4*a~2 + 4*a*b - 3*a + 2*b - 2 
if r>0: 

u=adjunta(s ,r ,u) 

return (u) 



t32=l; 
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t43=l; 

for p in primes (2, 500000) : 
if p°/.3==2: 

t32=t32*p 
if p%4==3: 
t43=t43*p 



def cri (q) : 

j=-l 
h=0 

while h==0 and j<wl-l: 
j=j+l 

if q°/.(w[j] [0])==w[j] [1] : 
h=l 
return (h) 

def f asf (q) : 
a=0 
h=l 

if cri(q)==0: 

if is_square(4*q+5)==0: 
if gcd(2*q+3,t43)==l : 
if gcd(q+2,t32)==l: 

if gcd(3*q+4,t32)==l: 
q0=q+l 
al=-l 
a0=-l 
h=0 

while h==0 and a0<(3*q/2) 
q0=q0+l 
al=al+4 
a0=a0+3 

for p in divisors (qO) 
if p°/„al==a0: 
h=l 

a=(a0-2)/3 

return ( [h,a] ) 
def cricu(q) : 
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if gcd(q+2,t3)==l : 

if gcd(3*q+4,t3)==l : 

if gcd(2*q+3,t4)==l : 
if gcd(q+3,t7)==l: 

if gcd(7*q+9,t7)==l : 
if (q+4)7„19>0: 

if (ll*q+14)7.19>0: 
if (q+5)7.11>0: 

if (15*q+19)7.11>0: 
if (q+7)7.17>0: 

if (23*q+29)7.17>0: 
h=0 
j=-l 

while h==0 and j<len(w)-l: 
j=j+l 

if q7.w[j] [0]==w[j] [1] : 
h=l 

if h==0: 

v. append (q) 

return (v) 



t=2*3*5*7*l 1*13* 17*19 
t 3=2*5*ll*17 
t4=3*7*ll*19 
t7=5*19 

w=cribata(t,2000) 
wl=len(w) 
print wl 



v=[] 

for q in xsrange(t): 

cricu(q) 
vl=len(v) 
print vl 
save(v,"qfl9") 

t=2*3*5*7*l 1*13* 17*19 
tc=t*23 

t3=2*5*ll*17*23 

t4=3*7*ll*19*23 

w=cribata(tc,2000) 

wl=len(w) 

print wl 

v7=v 

v=[] 

v2=tc/t 
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aO=-l 

for 11 in xsrange(0,v2) : 

for i in xsrange(0,vl) : 

q=v7[i]+ll*t 

cricu(q) 
vl=len(v) 
print vl 
save(v,"qf23") 

a0=0 

for 11 in xsrange (0,9000) : 
for i in xsrange (0, vl ) : 
ql=v[i]+ll*tc 
wt=fasf (floor (ql)) 
a=wt [1] 
if wt [0] ==0 : 

print ql," falla " 

a=0 
if a>a0: 

a0=a 

print ql,a 

6.2. A computer program verifying the q— Conjecture and Erdos-Straus 
Conjecture of Type-I. According to Lemma [5] if for each q G Z>o there exist 
x, y, z in Z>o such that one of the following relations holds 

q = 1 + 3x + 3y + Axy, 

q = 5 + 5x + 5y + Axy, 

q = {{Ay + 3)z + (3y + 2)){x + 1) - (y + 2). 

the Erdos-Straus Conjecture is true. We conjecture that all the natural number 
q G Z>o can be written as one of the three above relations. We refer it as the 
q— Conjecture. Observe that this is equivalent to for all prime numbers of the form 
4g + 5 the corresponding value q can be write as q = {{Ay + 3)z + {3y + 2))(x + 1) — 
{y + 2) for some x, y, z in Z>o- It is worst to observe that the identity 

4 1 1 1 

abc-a-b ~ a^- + a{ac-l)^ + (ac-l)^n 

with n — abc — a — 6, implies that — has a decomposition as sum of three Egyptian 
fraction with exactly one of the denominators multiplies of n (a decomposition of 
type I). We are verified q— Conjecture for q < 2.5 x 10 13 . For this object, we do 
a litter change in the program above. We write two new functions crice and fase 
instead of cricu and fasf, respectively. 

The Sage code of these function is: 
def crice (q) : 

if gcd(4*q+5,tc)==l : 
if gcd(q+2,t3)==l : 

if gcd(3*q+4,t3)==l : 

if gcd(2*q+3,t4)==l : 



18 



M. BELLO-HERNANDEZ, M. BENITO, AND E. FERNANDEZ 



if gcd(q+3,t7)==l : 

if gcd(7*q+9,t7)==l : 
if (q+4)7„19>0: 

if (ll*q+14)%19>0: 
if (q+5)7.11>0: 

if (15*q+19)7.11>0: 
if (q+7)7.17>0: 

if (23*q+29)7.17>0: 
h=0 
j=-l 

while h==0 and j<len(w)-l: 
j=j+l 

if q7.w[j] [0]==w[j] [1] : 
h=l 

if h==0: 

v. append (q) 



return (v) 



def f ase (q) : 
a=0 
h=l 

if cri(q)==0: 

if gcd(4*q+5,tf )==1: 

if gcd(2*q+3,t43)==l : 
if gcd(q+2,t32)==l: 

if gcd(3*q+4,t32)==l: 

if is_prime (4*q+5)==l : 
qO=q+l 
al=-l 
aO=-l 
h=0 

while h==0 and a0<(3*q/2): 
qO=qO+l 
al=al+4 
a0=a0+3 

for p in divisors (qO) : 
if p%al==aO: 
h=l 

a=(a0-2)/3 

return ( [h,a] ) 

6.3. A computer program checking the Erdos-Straus Conjecture. We are 

verified that all the natural numbers n, 2 < n < 2 x 10 14 , satisfy the Edos-Straus 
Conjecture. To this end we do a more sharp sieve than in the above programs. The 
new cricue function replace crice function and add several new testing taking into 
account Lemma [9] which gives a characterization for the convergence of the greedy- 
type algorithm. Some of these new checkup are given by the minimal number of 
steps needed for the convergence of the greedy-type algorithm. Moreover, we divide 
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the interval [2.5 x 10 13 , 5 x 10 13 ] in several intervals each one of 10 3 numbers, then 
using the function ferti, which uses fase function, we check the ESC in each interval. 
Next we show the Pari GP code of our program: 

f ase(q,h,q0,al,a0,p,q7)={ 
h=0; 

if (gcd(4*q+5,tp)>l, 
h=l, 

if (gcd (2*q+3 , t43) >1 , 
h=l, 

if (gcd(q+2,t32)>l, 
h=l, 

if (gcd(3*q+4,t32)>l, 
h=l, 

if (cri(q)==l, 
h=l; 
a0=2, 

if (ispr ime (4*q+5) ==0 , 
h=l; 
a0=2, 

q7=(3*q)\2; 

qO=q+l; 

al=-l; 

aO=-l; 

h=0; 

while (h==0&&a0<q7 , 
q0++; 
al+=4; 
a0+=3; 

f ordiv(qO,p, 

if (p°/.al==aO , 

h=l, ))))))))); 

[h,(a0-2)\3]} 

adjimta(s,r,u,c,h)={h=0;f or (i=l ,matsize(u) [1] ,if (s°/ u [i , 1] ==0 , 
c=gcd(s,u[i,l]) ;if (r°/.c==u[i,2] , 
h=l ,),)); if (h==0,u=concat (u, [s,r] ) , ) ;u} 

{u=matrix (3 , 2 , x , y , 0) ; u [1 , 1] =2 ; u [2 , 1] =3 ; u [3 , 1] =7} 
cr ibat a (t , m , b , a , s , r ) ={ 
f or(b=0,m, 
f or(a=0,b, 

s=4*a*b+3*a+3*b+2 ; 
if (t°/„s==0, 

r=4*a*b+2*a+3*b;u=adjunta(s ,r ,u) ,) ; 
s=4*a*b+3*a+3*b+2 ; 
if (t°/„s==0, 

r=4*a*b+3*a+2*b;u=adjimta(s,r ,u) ,) ; 
s=4*a*b+3*a+4*b+2 ; 
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if (t%S==0, 

r=3*a*b+2*a+3*b; 

u=adjunta(s,r,u) ,) ; 
s=4*a*b+4*a+3*b+2 ; 
if (t%s==0, 

r=3*a*b+3*a+2*b; 

u=adjunta(s,r,u) , ) ; 
s =4*a*b+3*a+4*b+3 ; 
if (t7 s==0, 

r=3*a*b+2*a+2*b ; 

u=adjunta(s,r,u) , ) ; 
s =4* a *b+4*a+3*b+3 ; 

if (ty oS ==o, 

r=3*a*b+2*a+2*b; 

u=adjunta(s,r,u) , ) ; 
s =4* a *b-4*a~2+3*b+2 ; 
if (t7„s==0, 

r=4*a~2*b-4*a~2-4*b~2*a+4*a*b-a-3*b~2+b; 

if (r>0, 

u=adjunta(s,r,u) ,) ,) ; 
s =4*a*b-4*a~2+3*b+2 ; 
if (t7„s==0, 

r=4*a~2*b-4*a~2-4*b~2*a+4*a*b+a-3*b~2; 

if (r>0, 

u=adjunta(s,r,u) ,) ,) ; 
s=4*a*b-4*a~2+4*b+2-a; 
if (t°/,s==0, 

r=4*a~2*b-3*a~2-4*b~2*a+4*a*b-a-4*b~2+b; 

if (r>0, 

u=adjunta(s,r,u) ,) ,) ; 
s =4*a*b-4*a~2+a+3*b+3 ; 
if (t%s==0, 

r=4*a~2*b-3*a~2-4*b~2*a+2*a*b-3*b~2-b; 

if (r>0, 

u=adjunta(s,r,u) ,) ,) ; 
s=4*a*b-4*a~2+a+3*b+2 ; 
if (t%s==0, 

r=4*a~2*b-3*a~2-4*b~2*a+2*a*b-3*b~2+a; 

if (r>0, 

u=adjunta(s,r,u) ,) ,) ; 
s =4*a*b-4*a~2-a+4*b+3 ; 
if (t%s==0, 

r=4*a~2*b-3*a~2-4*b~2*a+4*a*b-4*b~2-b; 

if (r>0, 

u=adjunta(s,r,u) ,) ,) ; 
s=-4*a~2+4*a*b-4*a+3*b-l ; 
if (s>0, 

if (t7.s==0, 
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r=- 4*a~2 + 4*a*b - 5*a + 3*b - 3; 
if (r>0, 

u=adjunta(s,r,u) ,),),); 
s=-4*a~2+4*a*b-4*a+3*b-l ; 
if (s>0, 

if (t7„s==0, 

r= - 4*a~2 + 4*a*b - 3*a + 2*b - 2; 
if (r>0, 

u=adjunta(s,r,u) ,),),))) ;u} 



cricue(q,h, j)={ 
h=l; 

if (gcd(4*q+5,tc)==l, 

if (gcd(q+2,t3)==l, 

if (gcd(3*q+4, t3)==l , 

if (gcd(2*q+3,t4)==l, 

if (gcd(q+3,t7)==l, 

if (gcd(7*q+9,t7)==l, 

if ((q+4)7„19>0, 

if ((ll*q+14)y„19>0, 

if ((q+5)7„ll>0, 

if ((15*q+19)7,ll>0, 

if ((q+7)°/„17>0, 

if ((23*q+29)7„17>0, 

h=0; 

j=0; 

while (h==0 && j<matsize(w) [1] , 

if ((q7.w[j,l])==w[j,2] , 
h=l,)),),),),),),),),),),),),);hj 

{t32=l ; t43=l ; f orprime (p=2 , 500000 , if (p7.3==2 , 

t32*=p,);if (p%4==3, 

t43*=p,))} 

cri(q, j ,h)={ 
j=0; 
h=0; 

while (h==0 && j<wl, 

if (q'/.w[j,l]==w[j,2] , 
h=l,));h> 



{tc=2*3*5*7*l 1*13*17* 19; 
t=tc ; 
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t3=2*5*ll*17; 
t4=3*7*ll*19; 
t7=5*19; 

w=cribata(t,5000) ; 
wl=matsize (w) [1] ; 
print (wl)} 
{v=[] ; 

for(l=0,t\6-l, 
q=6*l+5; 
h=cricue(q) ; 
if (h==0, 
k=0; 

while (h==0&&k<5, 
k++; 

r=((4*q+5)* (q+l+k) )%(4*k-l) ; 
c=((4*q+5)*(q+l+k))\(4*k-l) ; 
if (r==0, 
h=l, 

if ( ( (4*q+5) * (q+l+k) * (c+1) )7„ (4*k-l-r) ==0 , 
h=l,))),); 

if (h==0 , v=concat (v, q) , ) ) ; 
vl=matsize(v) [2] ; 
print ( " paso I ", vl)} 
{vl9=v; 

t=2*3*5*7*l 1*13*17*19; 

tc=t*23*9; 

t3=2*5*ll*17*23; 

t4=3*7*ll*19*23; 

w=cribata(tc , 5000) ; 

wl=matsize(w) [1] ; 

print (wl) ; 

x=0 ; v23=vector (4000000 , x) ; j 23=0 ; 
v2=tc/t;> 

{for(ll=0,v2-l, 

print 11; 
for(il9=l,vl, 
q=vl9[il9]+ll*t; 
h=cricue(q) ; 
if (h==0, 
k=5; 

while (h==0 && k<7, 
k++; 

r=((4*q+5)*(q+l+k))7„(4*k-l) ; 
c=((4*q+5)* (q+l+k) -r)\(4*k-l) ; 
if (r==0, 
h=l, 

if ( ( (4*q+5) * (q+l+k) * (c+1) ) % (4*k-l-r) ==0 , 
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h=l, ))),); 

if (h==0, 

j23++;v23[j23] =q,)))j 
x=0;v=vector(j23,x) ; 
for(i23=l, j23,v[i23] =v23[i23]) ; 
vl=matsize (v) [2] ; 
print("erti23 ",vl); writeCv.txt" ,v)> 

{tp=l ; t32=l ; t43=l ; f orprime (p=2 , 15000 , tp*=p ; 
if (p'/.3==2 , t32*=p , ) ; if (p°/„4==3 , t43*=p , ) ) } 

{t=2*27*5*7*l 1*13*17* 19*23; 
w=cribata(t , 100) ; 
wl=matsize (w) [1] ; 
print (wl)} 



ferti(li,lf ,a0,10,l,i7,q,y,a)={a0=0;10=matsize(v) [2] ; 
for(l=li,lf ,print(l) ;for(i7=l,10,q=v[i7]+t*l;y=fase(q) ;a=y[2] ; 
if (y[l]==0 .print (q," falla") ; write ( "salida23" ,q, "f alia") , 
if (a>a0,a0=a; print (" ",q," ",a); 
write ("salida23"," ",q," " , a) ,))))} 
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